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ABSTRACT 



Context. The Small Magellanic Cloud (SMC) is ideally suited to investigating the recent star formation history from X-ray source 
population studies. It harbours a large number of Be/X-ray binaries (Be stars with an accreting neutron star as companion), and the 
supernova remnants can be easily resolved with imaging X-ray instruments. 

Aims. We search for new supernova remnants in the SMC and in particular for composite remnants with a central X-ray source. 
Methods. We study the morphology of newly found candidate supernova remnants using radio, optical and X-ray images and inves- 
tigate their X-ray spectra. 

Results. Here we report on the discovery of the new supernova remnant around the recently discovered Be/X-ray binary pulsar 
CXOJ012745. 97-733256. 5 = SXP 1062 in radio and X-ray images. The Be/X-ray binary system is found near the centre of the 
supernova remnant, which is located at the outer edge of the eastern wing of the SMC. The remnant is oxygen-rich, indicating that 
it developed from a type lb event. From XMM-Newton observations we find that the neutron star with a spin period of 1062 s (the 
second longest known in the SMC) shows a very high average spin-down rate of 0.26 s per day over the observing period of 18 days. 
Conclusions. From the currently accepted models, our estimated age of around 10000-25000 years for the supernova remnant is not 
long enough to spin down the neutron star from a few 10 ms to its current value. Assuming an upper limit of 25000 years for the age 
of the neutron star and the extreme case that the neutron star was spun down by the accretion torque that we have measured during 
the XMM-Newton observations since its birth, a lower limit of 0.5 s for the birth spin period is inferred. For more realistic, smaller 
long-term average accretion torques our results suggest that the neutron star was bom with a correspondingly longer spin period. This 
implies that neutron stars in Be/X-ray binaries with long spin periods can be much younger than currently anticipated. 

Key words, galaxies: individual: Small Magellanic Cloud - stars: neutron - stars: Be - X-rays: binaries 



1. Introduction 

One of the most cataclysmic events in the universe is the ex- 
plosion of a massive star as a supernova, which can create a 
neutron star (NS). NSs are thought to be born rapidly spin- 
ning with rotation periods of a few 10 ms. Their rotation is 
first slowed down by magnetic dipole braking and then by the 
propeller effect when the NS is in a binary star system and 
mass loss from the companion star beg ins (for a review see e.g. 
iBhattacharva & van den Heuvell Il99l l When rotating slowly 
enough, accretion onto the NS sets in (Ghosh & Lamb 1978) and 
the system can be detected as an X-ray binary. The spin evolution 
of the NS in a high-mass X-ray binary (HMXB) system depends 
on the initial magnetic field stre ngth of the NS and on the mass 
accretion rate (lUrpin et al.'l998). 

Be/X-ray binary (BeXRB) systems are a subgroup of 
HMXBs with an NS accr eting matter from the circumstellar disk 
of a Be star ('Reig"201l|). The Small Magellanic Cloud (SMC) 
harbours an extraordina rily large number of BeX RB s (C oe et alj 
12005': 'Shtvkovskiv & Gilfanov 2005' 'Hab erl et alJl2008h . These 
are believed to have been creat ed during a burst of s tar forma- 
tion about 42 million years ago dAntoniou et alj|20Tol) . For more 



than fifty BeXRBs in the SMC the NS spin period is known 
from the detection of coherent pulsations in their X-ray flux. 
The spin periods ra nge from 2.16 s (XTE JOl 19-73 1 = SXP 2.16; 



Corbet etal.ll2003h to 1320 s (RX JO 103. 6-7201 = SXP 1323; 



Haberl & Pietsch 2005 ) with a bimo dal structure indicated in 
their distribution (Knigge et alj|201 lb . 



From ROSAT HRI data, iHughes & SmithI (11994 proposed 
two BeXRBs within supernova remnants (SNRs) in the SMC. 
For one of them, loca ted in the direction of the SNR IKT21 
(Ivan der Hevden et al. 2004), pulsations of 345 s were discov- 



* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by ESA 
Member states and the USA (NASA) 



ered bv llsrael et al.1 (l2000r SAX J0103.2-7209). The other - still 
unconfirmed - candidate BeXRB in the direction of IKT 25 was 
not detected by XMM-Newton or Chandra. In both cases it is 
not clear whether the BeXRB (or candidate) is physically asso- 
ciated with the SNR because they are not located near the centre 
of the remnant and the high number of BeXRBs in the bar of the 
SMC make chance coincidences likely. The SNR, which was 
created by the supernova explosion, is expected to have faded 
beyond detectability before accretion onto the NS starts. Finding 
a BeXRB associated to a SNR would therefore suggest a much 
faster spin-down than predicted by the standard models, or else 
the NS star was born with slow rotation. Here, we present the dis- 
covery of an SNR in radio and X-rays around the BeXRB pulsar 
CXOJ012745.97 -733256.5 = SXP 1062. The discovery of the 
new SMC pulsar and its possible association with a shell nebula 
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Fig. 1. Images of the region around the BeXRB SXP 1062. Left: Continuum-subtracted MCELS images. Red, green, and blue 
correspond to Ha, [Om] and [Sn]. The cross marks the centre of the white circle (R.A. = 01*'27"44';i5, Dec. = -73°33'01'.'6, 
J2000), which indicates the estimated size of the SNR in all images. The Be star counterpart of SXP 1062 is visible close to the 
cross. Middle: MOST 36 cm radio image, logarithmically scaled from 0.5 to 3 mJy. Right: Combined, instrument background- 
subtracted EPIC-pn and MOS colour image obtained from the XMM-Newton observations. Red, green, and blue denote X-ray 
intensities in the 0.2-1.0, 1.0-2.0, and 2.0-4.5 keV bands. The images from the individual bands were adaptively smoothed with a 
Gaussian cr of 4" (at high intensities) to 20" (low intensities). Circles indicate the extraction regions for the X-ray spectra: Source 
region for SXP 1062 (A, 29" radius) and associated background (B, 25"). For the SNR spectrum (C, 84.5"), regions D (40"), E 
(11"), and F (14") were excluded and region G (60 ") was used to estimate the background. Regions E and F were excluded because 
they contain possible point sources in the direction of the SNR as found in source detection runs on the EPIC images. 



around it has re cently been reported by iHenault-Brunet et al] 
(1201 IL hereafter lHBllh . 



2. Multi-wavelength observations and data analysis 

2.1. Optical 

The Magellanic Cloud Emission Line SurvejQ was carried out 
from the 0.6 m University of Michigan/CTIO Curtis Schmidt 
telescope, equipped with an SITE 2048x2048 CCD. This re- 
sulted in a field of 1.35° x 1.35° with a pixel scale of 2.4" x 
2.4". The Magellanic Clouds were mapped in narrow wave- 
length bands corresponding to Ha, [Oiii], and [Sii], together 
with matching red and green continuum bands. All the images 
were continuum-subtracted, flux-calibrated, and assembled into 
mosaic images. A composite image in the three bands of the 
region around SXP 1062 is presented in Fig.lT] revealing a shell- 
ty pe nebu la dominated by the strong [O iii] line as also discussed 
bv lHBllL 

2.2. Radio-continuum 

The field of SXP 1062 was observed with the Australia 
Telescope Compact Amy (ATCA) on 2009-01-05 with amy 
configuration 6C, and on 2009-02-05 in array configuration 
EW352 (ATCA project C1869). The field was observed at wave- 
lengths of 20 cm and 13 cm. Each session was carried out in 
"snap-shot" mode, totalling ~3.5 hr of integration time over a 
12 hr period. The se data were combined with the relevant parts 
of the data from Wong et al.l (12011) and FiUpovic et al.1 (Il995h 
to improve u-v plane coverage. The 13 cm data proved unusable 
due to radio frequency interference. The 20 cm data was useful 
only when long baselines were excluded thus the final image has 



' MCELS: rhttp://www.ctio.noao.edu/~mceIs71 



a resolution of 145"x 131" with anrm.s. noise of 0.5 mJy/beam. 
We report an integrated flux density of 6 mJy at 20 cm. 

Using the image of Ye & Turtle (1993), we measured an in- 
tegrated flux density of 9mJy at 36 cm (see Fig.[T]i. Combining 
this measurement with the 20 cm measurement we derive a spec- 
tral index (5v v"^) of or = -0.8 + 0.4 which is steeper than 
most SNRs in the Magellanic Clouds (Filipovic et al.lll998i) . We 
estimate that the surface brightness of this SNR is 1.6 x 10"^^ 
Wm^^Hz"' Sr"' with a luminosity of 3.4 x 10^^ WHz"', as- 
suming a distance of 60 kpc and diameter of 166". This makes it 
flie faintest known radio SNR in the SMC (.Filipovic et al...2005l 
l2008l:IOwenetani2011l) . 

2.3. X-rays 

Chandra (between 2010-03-31 and 2010-04-29) and XMM- 
Newton (2010-03-25 to 2010-04-12) observed the field around 
the star-forming region NGC 602 and revealed the new BeXRB 
SXP 1062 in the eastern wing of the SMC (HBll) with a 
spin period of 1062 s, the second longest known in the 
SMC after SXP 1323. Here, we present the analysis of the 
four XMM-Newton observations (observation IDs 0602520401, 
0602520201, 0602520301, and 0602520501) with durations of 
66.3, 119.0, 91.2, and 55.0 ks and spread over 18 days. 

2.3.1. Morphology 

Following the analysis of the SMC survey with XMM-Newton 
(Haberl et al. in preparation, Sturm et al. in prepara tion), we pro- 
duced X-ray images fro m the combined EP IC-pn dStriider et al.l 
12001 1) and EPIC-MOS dTurner et al .11200 lb data in the energy 
bands 0.2-1.0 keV, 1.0-2.0 keV, and 2.0-4.5 keV. Our analy- 
sis of the EPIC data used the tools from the XMM-Newton SAS 
package version 1 1 .0. 10 To optimise the sensitivity for faint dif- 

^ Science Analysis Software (SAS), [http://xmm.esac.esa.int/sas/| 
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fuse X-ray emission we removed intervals of high background 
caused by soft protons and subtracted the contribution of the de- 
tector particle background from the images. For EPIC-pn also 
the out-of-time events recorded during CCD readout were sub- 
tracted. The final colour image of the region around SXP 1062 is 
presented in Fig.[T] Although the emission from SXP 1062 dom- 
inates the hard energy bands, soft extended emission is clearly 
seen around the X-ray pulsar The extent of this emission is 
similar to the radio and MCELS images. We also checked the 
Chandra images for the soft emission, but the sensitivity at low 
energies is insufficient for the Chandra ACIS-I instrument to al- 
low a significant detection. 

2.3.2. Spectral analysis 

We extracted EPIC spectra from the X-ray pulsar and the soft 
emission region, merging the data from all four observations. 
The extraction regions are indicated in Fig. [1] and described be- 
low. No MOS 1 data is available because the source was located 
on CCD 6, which is no longer operative. For the spectra we 
selected single and double-pixel events (PATTERN 0-4) for pn 
and single to quadruple events (PATTERN 0-12) for MOS with 
FLAG=0 to avoid bad pixels. Applying the same background 
flare screening as for the images resulted in net exposures for 
the spectra of 206.8 ks for pn and 236.7 ks for MOS2. 

To extract the spectra of SXP 1062 we used a circular re- 
gion and a very close background region inside the soft emis- 
sion region to remove the contribution of this component. We 
fitted the pn and MOS2 simultaneously with a power-law, at- 
tenuated by two absorption components. The first accounts for 
the Galactic foreground absorption with a fixed column den - 
sity of 6x10^° cm and solar abundances dWilms et aUlIOOOl) . 
while the second models the absorption along the line of sight 
within the SMC and local to the source with free column den- 
sity and reduced abundances o f 0.2 for elements heavier than 
helium (iRussell & Dopitalll992h . This resulted in an acceptable 
fit with reduced of 1.12 for 377 degrees of freedom (see 
Fig. |2]). The best-fit value for the column density in the SMC 
was 1.8 + 0.2 xlO^' cm"^ and the power-law p hoton index 
0.74 ± 0.02, values typical of BeXRBs in the SMC faaberl et al.1 
120081) . The observed average flux (0.2-10 keV) during the 
XMM-Newton observations was 1.4x10"'^ erg cm"^ s"' (as ob- 
tained from the pn spectrum), which corresponds to a source lu- 
minosity corrected for absorption of 6.3x10"'^ erg s"' (assum- 
ing a distance of 60 kpc). The flux increased after the first ob- 
servation and finally dropped by 27% in the last observation 
(1.33x10-12 erg cm-2 s"', 1.53x10-'^ ergcm^^ ^-i^ 1.54x10-'^ 
erg cm^^ s"' and 1.13x10"'^ erg cm"^ s"' in chronological or- 
der), which could indicate that we w itnessed a typical type I out- 
burst dOkazaki & Negueruelall2001h . 

To investigate the soft extended emission, we extracted the 
EPIC spectra from a ring around the position of SXP 1062. The 
relatively wide wings in the point spread function of the XMM- 
Newton telescopes means a contribution from the X-ray pulsar to 
the soft emission cannot be avoided. Therefore, we modelled the 
spectrum of the diffuse emission with a power-law component in 
addition to a thermal component. We kept the SMC absorption 
and photon index for the power-law fixed at the values derived 
from the spectra of SXP 1062, allowing only a free normalisa- 
tion. For the thermal component we used the plasma emission 
model mekal (also with SMC abundances as for the absorption) 
available in XSPEC (Meweetal. 1985) with SMC absorption 
(free in the fit and different to that of the power-law compo- 
nent). The resulting best-fit parameters (reduced ;t'2 = 1.22 for 
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Fig. 2. EPIC-pn (black) and MOS2 (red) spectra of SXP 1062 
combining the data from the four XMM-Newton observations. 
The best-fit ower-law model is plotted as histogram, and the 
lower panel shows the fit residuals. 
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Fig. 3. EPIC spectra as in Fig.|2]extracted from the soft emission 
region around SXP 1062. The best-fit model is composed of a 
thermal (dashed line) and a power-law (dotted) component (see 
text). 



16 degrees of freedom) are 7.3+^-^x102" cm for the SMC ab- 
sorption and 0.23+0.05 keV for the temperature. The observed 
flux (0.2-2 keV) of 1.1x10"''* erg cm"^ g-i makes it the X-ray 
faintest SNR known in the SMC (for other faint remnants see 
iFilipovic et alj|2008l) . The spectra with the best-fit model are 
presented in Fig. [3] Unfortunately, the low-quality statistics of 
the spectra do not allow a more detailed comparison of different 
models for the thermal spectral component. 

The Hi column density of the SMC i n the direction of 
SXP 1062 is measured to 2.1xl02' cm ^ dStanimirovic et all 
1999). The absorption inferred from the X-ray spectrum of 
SXP 1062 is similar, while for the soft emission region we can 
only place an upper limit of ^1.6x10^' cm 2. The fact that 
BeXRBs show source intrinsic absorption and the large uncer- 
tainty in the column density of the soft X-ray emission spectrum 
do not allow deriving any constraints on the radial position of the 
objects within the SMC. Nevertheless, the two column densities 
inferred from the X-ray spectra of the BeXRB and the diffuse 
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X-ray emission region (the SNR) are consistent with a location 
of the BeXRB within the SNR. 



2.3.3. Timing analysis of SXP 1062 

We performed a timing analysis of the EPIC light curves using 
the Bayesian approach (Gresorv & Loredo 1996) and a R ayleigh 
7} tes t for one harmonic (iHaberl & Zavlinii2002i : .Buccheri et al.l 
Il983h around the known periodicity of 1062 s. To investigate 
a possible evolution of the spin period we first analysed each 
XMM-Newton observation separately. The obtained spin peri- 
ods are plotted versus time in Fig. |4] and show a steady spin- 
down trend over the observing interval of about 18 days. A lin- 
ear fit to the spin period evolution (line in Fig. |4| results in a 
period change P of 2.8+|}^xl0-*' s s"' (0.24+°{;^ s day-'). This 
value was used as starting point for a phase-coherent analysis 
combining the EPIC-pn data of all four observations searching 
a grid of different values for the spin period P and a constant P. 
We used the 7} test with one harmonic to account for the first 
harmonic seen in the FFT power density spectrum. This resulted 
in a maximum Zj value of 344.1 for P - 1061.24 s (at an epoch 
of MJD 55280.53, the start of the first EPIC-pn exposure) and 
P=3.0±0.5xlO-''ss-'. 



3. Discussion and conclusions 

Our multi-wavelength morphological studies of the field around 
the new BeXRB pulsar SX P 1062 confirm the previously 
unknown SNR proposed by IHB 1 ll Its SNR nature is fur- 
ther supported by its X-ray spectrum obtained from both the 
XMM-Newton EPIC data and its steep radio spectrum. The 
stron g [Oml fine emission sugg ests the SNR as oxygen-rich 
type (lArbutina & Urosevicll2005l) . which usually comes from a 
type lb event, the explosion of a massive O, B, or WR star. 

SXP 1062 is detected close to the centre of the SNR. At 
the outer edge of the eastern wing of the SMC, only very few 
BeXRB s are dete cted, and all known SNRs are located in the 
bar of the SMC ( Payne et allboOTt iFihpovic et a"l]|2008l) . This 
makes a chance coincidence of SXP 1062 with the SNR very 
unlikely. However, the ages of SNRs and BeXRBs are expected 
to be rather different. While SNRs can usually be seen for at 
most a few 10^ years in X-ray s and radio (for the methods used 
to estimate SNR ages see e.g. Xu et al. 2005), BeXRBs are ex- 
pected to become X-ray active much later (see below). Using the 
relation ty = 3.8x10^ Rpc(kT)^J^^ from IXuet all ( l2005h of tem- 
perature derived from X-ray spectral modelling (0.23 keV) and 
the size of the SNR (~2.5' in diameter, corresponding to ~40 
pc at a distance of 60 kpc), we estimate the age of the SNR to 
~16 ky. If on e compa res the ages e stimated for SNRs in the 
SMC by van d er Hevde iiet al.l (l2004t) with their sizes as listed in 
iBadenes et al.l (1201 0). one finds ages of at most 25 ky for rem- 
nants with diameters around 40-50 pc. 

We estimate the distance between the BeXRB and the cen- 
tre of the SNR to about 9", which coiTesponds to ~8xlO'^ km 
at SMC distance. From evolutionary investigations, kick veloc- 
ities around 100- 200 km s~' were inv oked for neutron stars in 
HMXB systems dPortegies Zwartll995h . Using this as maximum 
space velocities for the binary systems and neglecting a veloc- 
ity component in radial direction, it would take 13-26 ky for 
SXP 1062 to move from the SNR centre to its current position, 
compatible with the SNR a ge estim ates from above, which are 
also consistent with those of IHB 1 iL 
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Fig. 4. Evolution of the spin period of SXP 1062 as obtained 
from the combined EPIC-pn and EPIC-MOS2 data of the four 
XMM-Newton observations. Spin periods and 1-cr errors were 
derived with the Bayesian odds ratio method. 



Statistical measure ment s of Be XRBs in the SMC result in 
lower space velocities: ICod (l200l inferred a lower limit of 30 
km s ' from measuring the separation of BeXRBs to nearby 
young stellar clusters. Similarly Antoniou et al. (2010) argue 
from the correlation between the number of BeXRBs and the 
local star formation rate 42 My ago at their current position for a 
maximum velocity of 15-20 km s"' . Both cases should actually 
be regarded as minimum velocities as the BeXRB mi ght not have 
been b orn in the nearest star cluster The work of Knigge et al.l 
(|20TT1) might also suggest that BeXRB pulsars with long spin 
period received stronger kicks. In either case, using lower space 
velocities (and also components in radial direction) would cor- 
respondingly increase the SNR age estimate^ 

According to the model calculations of lUrpin et al.l (Il998l) . 
it seems very unlikely that an NS can be spun down to a period 
around 1000 s within a few 10"^ years if it was born with a spin 
period around 10 ms. Of course, such model calculations assume 
simplified approximations of the conditions in a BeXRB (in par- 
ticular concerning the highly variable accretion usually observed 
in such systems), but only for models with extreme conditions 
(high magnetic field strength for the NS and low accretion rates), 
spin periods in excess of a few 100 s are reached after more than 
lO*" years. If the NS in SXP 1062 is indeed the compact rem- 
nant from the supernova explosion, which also created the SNR 
around it, then either the NS was born with a much longer spin 
period than 10 ms, or some other effect plays an important role 
to efficiently decelerate the rotation of the NS. 

SXP 1062 was discovered as a transient BeXRB pulsar dur- 
ing XMM-Newton and Chandra observations in 2010 when 
its X-ray luminosity was at the level of 6x10"^^ erg s"'. The 
BeXRB was not active du ring the ROSAT PSPC observations 
of the SMC (IHaberl etal.ll2000) . Based on the longest PSPC 
observation from 1991-10-07 and assuming the spectral shape 
measured by EPIC, we derive an upper limit of 1.5x10 '"^ erg 
cm"^ s"' (0.2-10 keV), assuming an upper limit of 10 cts and 
a vignetting-corrected exposure of 9.7 ks at the position of 
SXP 1062. This is about a factor of 10 lower than the flux dur- 
ing the XMM-Newton observations. While BeXRB transients 
generally exhibit spin-up during intervals of increased accretion 
(e.g. Bildsten et al. 1997), SXP 1062 is peculiar as it showed a 
strong spin-do wn consis tent with a constant rate of 0.26 s day ' 
over 18 days. ICoe et al ](|2010) investigated RXTE monitoring 
data covering ~10 years to study the spin period changes of 
BeXRB pulsars in the SMC over their active periods, which typ- 
ically last 50-500 days. They find four out of 15 systems with 
measured short-term spin period changes, which showed spin- 
down during outburst. In Table [1] we summarise the pulse pe- 
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Table 1. X-ray binary pulsars with periods of spin-down. 



Reference 

Bildsten et al. (1997) 
Cpe et al. ( 2010') 
_Makishima et al. (1988) 
Coe et al. (2010) 
this work 
Coe et al. (2010) 



Object 


Pulse 


f 


name 


period [s] 


[Hz s-'] 


Cen X-3 


4.82 


-(l-5)xlO-i2 


SXP 8.80 


8.90 


-5.1x10^'' 


GX 1+4 


110.2 


-2.1x10-" 


SXP 144 


144.5 


-5.9xlO-'3 


SXP 1062 


1062 


-2.6xl0-'2 


SXP 1323 


1325 


-5.0x10-'^ 



riod changes for some X-ray binary pulsars, which show times 
of spin-down. We list the frequency derivative f = -P/P^, which 
is directly related to the accretion torque exerted on the NS. We 
include the three SMC BeXRB pulsars with well measured val- 
ues of P and for comparison Cen X-3 (a super-giant HMXB 
with Roche-lobe overflow) and GX 1 +4, a low-mass system with 
M5 III donor star SXP 1062 exhibits spin-down about five times 
larger than the other listed SMC BeXRBs (which have remark- 
abl y similar values) . It should be noted that the values derived 
bv lCoe etal.ld2010h are average values over longer periods of 
time and that intermittent intervals of spin-up (as usually seen 
from accreting pulsars, Bildsten et al. 1997) reduce the longer- 
term average. 

SXP 1062 shows a remarkably large accretion torque, sim- 
ilar in magnitude to that of Cen X-3. In this HMXB intervals 
of steady spin -up and spin-down a lternate, which last typically 
10-100 days dBildsten et al.lll997h . but at much higher X-ray 
luminosities of ~10^^ erg s~' (e.g. iNagase et al.l Il992h . Also 
GX 1+4 was still observed at 10 ^^ erg s ' during its extended 
lowstate (Makishima et al. 1 988h . about a factor of 2 brighter 
than SXP 1062. This demonstrates that the generally accepted 
model for accretion, where the torque is directly proportional to 
the mass accretion rate (hence luminosity) is probably too sim- 
ple, as similar torques can be exerted, although the luminosity is 
more than a factor of 100 different. The work of Bildsten et al. 
(Il997l) suggests that disk-accreting pulsars are subject to in- 
stantaneous torques with similar magnitude but opposite sign. 
Different average long-term spin-up or spin-down values of in- 
dividual pulsars would then be the result of different time scales 
for reversals between spin-up and -down. 

SXP 1062 most likely also shows intervals of alternating 
spin-up and spin-down, which would reduce the magnitude of 
the long-term spin change. Assuming the extreme case that the 
NS was spun down with -2.6x10"'^ Hz s"' (-8.4x10^^ Hz y"') 
over its whole life of 25 ky (the maximum age) an upper limit of 
2 Hz for the total spin-down is derived, which corresponds to a 
lower limit of 0.5 s for the spin period at birth. For a more real- 
istic, smaller long-term average accretion torque a coiTespond- 
ingly longer birth spin period is expected. Therefore, if the NS 
in SXP 1062 is indeed the compact remnant from the supernova 
explosion that created the SNR, our results show that the NS was 
most likely born with a spin period much longer than a few tens 
of ms as generally adopted for an NS at birth. 

We confirm the existence of a new SNR around the 
BeXRB SXP 1062 in the SMC, independently suggested by 
iHenault-Brunet et al.l (1201 lb . The SNR is detected in optical, ra- 
dio, and X-ray wavelengths and shows a shell-like structure with 
the BeXRB close to its projected centre. We estimate an age of 
around 10-25 ky for the remnant. Such a time scale is too short 
to spin down the neutron star in a BeXRB to ~ 1000 s if it was 



born with a canonical spin period of a few 10 ms. SXP 1062 
is remarkable for showing a strong average spin-down rate f of 
-2.6x10 Hz s"' during the XMM-Newton observations dis- 
tributed over 18 days, while it was detected during high X-ray 
activity at 6.3x10-^^ erg s"' . Assuming the extreme case that such 
a high reverse accretion torque was exerted on the neutron star 
for its whole life, a lower limit of 0.5 s for the birth spin period 
is inferred. If the neutron star in SXP 1062 is indeed the compact 
remnant from the supernova explosion that created the SNR, we 
conclude that the neutron star in SXP 1062 was born with an 
even longer spin period for more realistic, lower long-term aver- 
age accretion torques. 
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